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Abstract 

We study the singular limit of a system of partial differential equa- 
tions which is a model for an aggregation of amoebae subjected to three 
effects: diffusion, growth and chemotaxis. The limit problem involves 
motion by mean curvature together with a nonlocal drift term. We 
consider rather general initial data. We prove a generation of interface 
property and study the motion of interface. We also obtain an optimal 
estimate of the thickness and the location of the transition layer that 
develops 

1 Introduction 

Let us start by a short description of life-cycles of the cellular slime molds 
(amoebae). The cells feed and divide until exhaustion of food supply. Then, 
the amoebae aggregate to form a multicellular assembly called a slug. It 
migrates to a new location, then forms into a fruiting body, consisting of 
a stalk formed from dead amoebae and spores on the top (fruiting bodies 
that are visible to the naked eye are often referred to as mushrooms) . Under 
suitable conditions of moisture, temperature, spores release new amoebae. 
The cycle then repeats itself. 

It is known that the aggregation stage is mediated by chemotaxis, i.e. 
the tendency of biological individuals to direct their movements according 
to certain chemicals in their environment. The chemotactant {acrasin) is 
produced by the amoebae themselves and degraded by an extracellular en- 
zyme (acrasinase) . For more details on the biological background, we refer 
to dS], [20] or [TO]. 

So the amoebae have a random motion analogous to diffusion coupled 
with an oriented chemotactic motion in the direction of a positive gradient 
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of acrasin. In 1970, Keller and Segel [15] proposed the following system as 
a model to describe such movements leading to slime mold aggregation: 



(KS) 



Ut = V • {D2VU) - V • (DiVv) 
vt = D^Av + f{v)u — k{v)v, 



inside a closed region 0. Here, u, respectively v, denotes the concentration 
of amoebae, respectively of acrasin; f{v) is the production rate of acrasin, 
and k{v) the degradation rate of acrasin (due to acrasinase); D2 = D2{u, v), 
respectively Di = Di{u, v), measures the vigor of the random motion of the 
amoebae, respectively the strength of the influence of the acrasin gradient 
on the flow of amoebae; Dy is a positive and constant diffusion coefficient. 
The problem is completed by initial data uq and vq and, assuming that 
there is now flow of the amoebae or the acrasin across the boundary dil., by 
homogeneous Neumann boundary conditions 



being the unit outward normal to dil.. 
An often used simplified model is obtained as follows. By some receptor 
mechanism, cells do not measure the gradient of v but of some with a 
sensitive function x satisfying x' > 0, so that Di{u,v) = ux'{v). By taking 
D2, f and k as constant functions and using some rescaling arguments, the 
system reduces to 



with du, dy, T and 7 some positive constants. 

Many analyses of the Keller-Segel model for the aggregation process were 
proposed. Chemotaxis having some features of "negative diffusion" , Nanjun- 
diah [20] suggests that the whole population concentrates in a single point; 
we refer to this phenomenon as the chemotactic collapse. In mathematical 
terms, this amounts to blow up in finite time. As a matter of fact, it turns 
out that the possibility of collapse depends upon the space dimension. In 
particular it never happens in the one-dimensional case whereas in two space 
dimensions, assuming radially symmetric situations, it only occurs if the to- 
tal amoebae number is sufficiently large. The problem of global existence 
and blow up of solutions has been intensively studied; we refer in particular 



In a different framework, Mimura and Tsujikawa [T7j, consider aggregat- 
ing pattern-dynamics arising in the following chemotaxis model with growth: 



Vn • V = Vv ■ u = on do, X {0, +00) 




Ut = duAu - V • (nVx(v)) 
TVt = dyAv + u — jv, 



to [9], [22], m, m, m, m, m- 




e^Au - eV • {uVx(v)) + f{u) 
Av + u — 7u, 
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where e > is a small parameter. The function / is cubic, and 1 being 
its stable zeros, and satisfies / > 0. In this model, the population is 
subjected to three effects: diffusion, growth and chemotaxis. The diffusion 
rate and the chemotactic rate are both very small compared with the growth 
rate. They observe that, in a first stage, internal layers — which describe the 
boundaries of aggregating regions — develop; in a second stage, the motion 
of the aggregating regions — which can be described by that of internal 
layers — takes place. The balance of the three effects (diffusion, growth and 
chemotaxis) makes the aggregation mechanism possible. Taking the limit 
e — > 0, they formally derive the equation for the motion of the limit interface 
and study the stability of radially symmetric equilibrium solutions. 

The purpose of this paper is to extend some of the results obtained by 
Bonami, Hilhorst, Logak and Mimura [1] about the singular limit of a variant 
of system (MT^), where the second equation is elliptic (r = 0): 



ut = Au-V- (nVx(f )) + i^Mu) mQx (0, +oo), 

= Av + u — jv in U X (0, +oo), 
du dv 

— = — = on 9ri X (0, +oo), 
ov ov 

u{x, 0) = uq{x) in 



in 



where is a smooth bounded domain in M (N > 2), is the Euclidian 
unit normal vector exterior to 90. We assume that 7 is a positive constant 
and that the nonlinearity/^ is given by 

fe{u) = u{l - u){u - ^) + eau{l - u) 

2 (1.1) 

=: f{u) + eg{u), 

with a > 0. The role of the function g is to break the balance of the 
two stable zeros slightly. The sensitive function % is smooth and satisfies 
x'{v) > for u > 0. 

We also assume that the initial datum satisfies uq £ C^(r2) and uq > 0. 
Throughout the present paper, we fix a constant Co > 1 that satisfies 

ll'"o|lco(T!) + l|V"o|lco(Q) + \\^uo\\^o(Q) < Co. (1.2) 

Furthermore we define the "initial interface" Tq by 

To := {x G n, Uo{x) = 1/2}. 

We suppose that Tq is a 0"^+^ hypersurface without boundary, for a i9 E 
(0, 1), such that, n being the Euclidian unit normal vector exterior to Tq, 

ToCCO and Vno(x) • n(x) / ifxETo, (1.3) 
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uo > 1/2 in Mo < 1/2 in (1.4) 

where Q^^^ denotes the region enclosed by Tq and the region enclosed 
between and Tq. 

The existence of a unique smooth solution to Problem (P^) is proved in 
[1], Lemma 4.2: 

Lemma 1.1. There exists eq > such that, for alls E (0, eo)? there exists a 
unique solution {u^,v^) to Problem (P^) on 0, x [0, +oo), with < < Cq 
on Qt- 

To study the interfacial behavior associated with this model, it is useful 
to consider a formal asymptotic limit of Problem (-P^) as e — > 0. Then the 
limit solution u^{x,t) will be a step function taking the value 1 on one side 
of the interface, and on the other side. This sharp interface, which we 
will denote by Ff, obeys a law of motion, which can be obtained by formal 
analysis (see Section [2|): 

Vn = -{N - 1)k + + V2a onFi, 

on 



-Am° + 7u° = u° inOx(0,r], 

— = onaOx(0,rl, 
I. ou 

where Vn is the normal velocity of Fj in the exterior direction, k the mean 
curvature at each point of F^. We set Qt := 0,x [0, T] and for each t G [0, T], 
we define q[^^ as the region enclosed by the hypersurface Ft and Qf'^ as the 
region enclosed between and Fj. The step function is determined 
straightforwardly from Ft by 

n . s f 1 in Q^^^ , , , , 

^i^^t) = {^ . Jo) fortG[0,r]. (1.5) 
10 m ill 

By a contraction fixed-point argument in suitable Holder spaces, the well- 
posedness, locally in time, of the free boundary Problem (P^) is proved in 
[1], Theorem 2.1: 

Lemma 1.2. There exists a time T > such that (P^) has a unique solution 
(i;0,F) on [0,T], with 



F= U (r* X {t}) G C2+'^ 



2+1? 
2 . 



0<t<T 



andv^lr e C^+^' 2 



2+1? 
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Bonami, Hilhorst, Logak and Mimura [1] have proved a motion of in- 
terface property; more precisely, for some prepared initial data, they show 
that {u^,v^) converges to {u^,v^) as e ^ 0, on the interval (0, T). So the 
evolution of Tt determines the aggregating patterns of the individuals. Here 
we consider the case of arbitrary initial data. Our first main result. Theo- 
rem [L3l describes the profile of the solution after a very short initial period. 
It asserts that, given a virtually arbitrary initial datum uq, the solution 
quickly becomes close to 1 or 0, except in a small neighborhood of the initial 
interface Tq, creating a steep transition layer around Tq [generation of in- 
terface). The time needed to develop such a transition layer, which we will 
denote by t^, is of order lne|. The theorem then states that the solution 

remains close to the step function on the time interval (motion 
of interface). Moreover, as is clear from the estimates in the theorem, the 
"thickness" of the transition layer is of order e. 

Theorem 1.3 (Generation and motion of interface). Let r] G (0,1/4) be 
arbitrary and set 

fi = f'{l/2) = 1/4. 

Then there exist positive constants Eq and C such that, for all e G (0,eo), 
all <t <T, where = lne|, we have 

'[-7?,l+r/] if xeMceiTt) 

u'{x,t) G I [-r^,r]] if x e ^]f^AAc,(^^) (1.6) 

+ if x€n^l\Afcs(rt), 

where Mr{Tt) ■= {x G il.,dist{x,Tt) < r} denotes the r -neighborhood of Ft- 



Corollary 1.4 (Convergence). As e ^ 0, the solution {u^,v^) converges to 

Jo<t<T(^t 



{u^,v^) everywhere in Uo<t<T(^t'' °^ ^ {*})• 



The next theorem deals with the relation between the set Tf := {x € 
Q,u^{x,t) = 1/2} and the solution Fj of Problem (P^). 

Theorem 1.5 (Error estimate). There exists C > such that 

TtcAfceirt) forO<t<T. (1.7) 

Corollary 1.6 (Convergence of interface). There exists C > such that 

dn{rirt)<Ce forO<t<T, (1.8) 

where 

d')-({A,B) := maxjsup (i(a, B), supd{b, A)} 

denotes the Hausdorff distance between two compact sets A and B. Conse- 
quently, rf — > as £ ^ 0, uniformly in < t < T, in the sense of the 
Hausdorff distance. 
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As far as we know, the best thickness estimate in the hterature was of 
order e|lne| (see [5], [6]). We refer to a forthcoming article [U], respec- 
tively [T], in which an order e estimate is established for a Lotka-Volterra 
competition-diffusion system, respectively for the FitzHugh-Nagumo sys- 
tem. 

The organization of this paper is as follows. Section [2] is devoted to 
preliminaries: we recall the method of asymptotic expansions to derive the 
equation of the interface motion; we also recall a relaxed comparison prin- 
ciple used in [1] . In Section [3l we prove a generation of interface property. 
The corresponding sub- and super-solutions are constructed by modifying 
the solution of the ordinary differential equation ut = e~'^f{u), obtained 
by neglecting diffusion and chemotaxis. In Section HJ in order to study the 
motion of interface, we construct a pair of sub- and super-solutions that rely 
on a related one-dimensional stationary problem. Finally, in Section El by 
fitting the two pairs of sub- and super-solutions into each other, we prove 
Theorem 11.31 Theorem 11.51 and theirs corollaries. 



2 Some preliminaries 
2.1 Formal derivation 

A formal derivation of the equation of interface motion was given in [3] . Nev- 
ertheless we briefly present it in a slightly different way: we use arguments 
similar to those in |21j where the first two terms of the asymptotic expansion 
determine the interface equation. The observations we make here will help 
the rigorous analysis in later sections, in particular for the construction of 
sub- and super-solutions for the study of the motion of interface in Section 

H 

Let {u^,v^) be the solution of Problem (P^). We recall that Tf := {x € 
O, u'^{x, t) = 1/2} is the interface at time t and call F^ := IJj>o(Ff x {t}) the 
interface. Let F = Uo<t<T(-'^t ^ i^i) solution of the limit geometric 

motion problem and let d be the signed distance function to F defined by: 

dist(a;,Ft) for x G , 



i'^''^'^;'! (2.1) 



where dist(x,F() is the distance from x to the hypersurface Tt in 0. We 
remark that d = on F and that |Vd| = 1 in a neighborhood of F. We then 
define 

q(!)= U {ni'^x{t}), qP= U (i^fxW). 

0<t<T 0<t<T 

We assume that the solution has the expansions 

u%x,t) = {0 or 1} + eui{x,t) -\ (2.2) 
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away from the interface T (the outer expansion) and 

u%x,t) = Uo{x,t,^^) + eUi{x,t,^^) + ■ ■ ■ (2.3) 

near T (the inner expansion). Here, the functions Uk{x,t, z), k = 0, 1, • • • , 
are defined for x £ i}, t > 0, z £ M. The stretched space variable ^ := 
d{x,t)/£ gives exactly the right spatial scaling to describe the rapid transi- 
tion between the regions {u^ ~ 1} and {u^ ~ 0}. We use the normalization 
conditions 

Uoix,t,0) = 1/2, Uk{x,t,0)=0, 

for all k > 1. The matching conditions between the outer and the inner 
expansion are given by 

Uo{x,t,+oo) = 0, Uk{x,t,+oo) =0, ^^^^ 
Uo{x,t, —oo) = I, Uk{x,t, —oo) = 0, 

for all A; > 1. We also assume that the solution has the expansion 

v''{x,t) =voix,t)+£Vi{x,t) + --- (2.5) 

in J7 X (0,r). 

We now substitute the inner expansion p. 31) and the expansion ()2.5p 
into the parabolic equation of (-P^) and collect the terms. We omit the 
calculations and, using \Vd\ = 1 near Tt, the normalization and matching 
conditions, we deduce that Uo{x,t,z) = Uo{z) is the unique solution of the 
stationary problem 

f uo" + fm = 

\ C/o(-oo) = 1, Uo{0) = 1/2, C/o(+oo) = 0. 

This solution represents the first approximation of the profile of a transition 
layer around the interface observed in the stretched coordinates. Recalling 
that the nonlinearity is given by f{u) = u{l — u){u — 1/2), we have 

%(.) = i(l-tanh^)=-£^. (2.7) 

We claim that Uq has the following properties. 

Lemma 2.1. There exist positive constants C and A such that the following 
estimates hold. 

< Uo{z) < Ce^^\^\ for z > 0, 
< 1 - Uq{z) < Ce-^l^l for z<0. 

In addition, Uq is a strictly decreasing function and 

\Uo'{z)\ + \Uo"{z)\ < Ce-^l^l for z G M. 
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Next we collect the terms. Since Uq depends only on the variable 
z, we have VUqz = which, combined with the fact that |V(i| = 1 near Tt, 
yields 

Ui,, + f'mUi = Uo'{dt -Ad + Vd- Vx{v^)) - 9{Uo), (2.8) 

a linearized problem corresponding to ()2.6p . The solvability condition for the 
above equation, which can be seen as a variant of the Fredholm alternative, 
plays the key role for deriving the equation of interface motion. It is is given 
by 

' = 0, 



Uo'\z){dt -Ad + Vd- Vx(t;°))(x, t) - g{Uo{z))Uo'iz) 

for all {x,t) G Qt- By the definition of g in (II. we compute 

g{Uo{z))Uo'{z)dz = - g{u)du = -a/6, 
; Jo 

whereas the equality (|2.7|) yields 

/ Uo'\z)dz = -= / -r—-^du = 1/6^2. 



Combining the above expressions, we obtain 

^dt - Ad + Vd- Vx{v^)) {x, t) = -V2a. (2.9) 



Since Vfi (= d(x, t)) coincides with the outward normal unit vector to the 
hypersurface Fj, we have dt{x,t) = —Vn, where Vn is the normal velocity of 
the interface F^. It is also known that the mean curvature k of the interface 
is equal to Ad/{N — 1). Thus the above equation reads as 

Vn = -{N - 1)k + + V2a on Ff, (2.10) 

on 

that is the equation of interface motion in (P^)- Summarizing, under the 
assumption that the solution -u^ of Problem (P^) satisfies 

u ^ < 7n\ as e ^ (J, 

\0 in qP 

we have formally showed that the boundary F^ between Q.f'^ and moves 
according to the law (j2.10p . 

One can note that, using the equality ()2.7p . we clearly have \/2aC/o' + 
g{Uo) = so that, substituting (^^9^ into i^B) yields Ui = 0. 
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2.2 A comparison principle 

The definition of sub- and super-solutions is the one proposed in [1]. 

Definition 2.2. Let {u~,uf) be two smooth functions with < uf on 
Qt and 

5ur dut „^ , 

< < — ^ on dn X (o,r . 

By definition, (uj,uf) is a pair of sub- and super-solutions if, for any 
which satisfies 



tig < —Av'^ + jv'^ < uf on Qt, 

— = ondnx(0,T), 
ov 



[2.11] 



we have 

where the operator L^je is defined by 

As proved in [1], the following comparison principle holds. 

Proposition 2.3. Let a pair of sub- and super- solutions be given. Assume 
that, for all x G 0, 

uj{x,0) < uo{x) < uf{x,0). 

Then, if we denote by {u^,v^) the solution of Problem (P^), the function 
satisfies, for all {x,t) G Qt, 

u~{x,t) < u'^{x,t) < uf{x,t). 

3 Generation of interface 

In this section we study the rapid formation of internal layers in a neigh- 
borhood of Tq = {x € ^1,uq{x) = 1/2} within a very short time interval of 
order lne|. In the sequel, we shall always assume that < t] < 1/4. The 
main result of this section is the following. 

Theorem 3.1. Let ij be arbitrary and define fi as the derivative of f{u) at 
the unstable equilibrium u = 1/2, that is 

= f'{l/2) = 1/4. (3.1) 

Then there exist positive constants eq and Mq such that, for all e G (0,eo), 
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• for all X € 

-r] <u''{x,n-'^e^\liie\) <l + r], (3.2) 

• for all X ^ Q such that \uq{x) — ^| > Mqe, we have that 

if uo{x)> 1/2 + Mqe then /^-^e^] lne|) > 1 - 77, (3.3) 

if uo{x) <l/2- Mqe then w^(a;, /x'^e^l lne|) < 77. (3.4) 

The above theorem will be proved by constructing a suitable pair of sub 
and super-solutions. 

3.1 The perturbed bistable ordinary differential equation 

We first consider a slightly perturbed nonlinearity: 

fs{u) = f{u) + 6, 

where 5 is any constant. For \6\ small enough, this function is still cubic and 
bistable; more precisely, we claim that it has the following properties. 

Lemma 3.2. Let 5q> Q he small enough. Then, for all 5 G {—5q,5q), 

• fs has exactly three zeros, namely a_((5) < a{6) < a+{5). More pre- 
cisely, 

fs{u) = {u- Q_(5))(Q+(5) - u){u - a{6)), (3.5) 
and there exists a positive constant C such that 

\a-{6)\ + \a{6) - 1/2| + \a+{6) - 1| < C\d\. (3.6) 

• We have that 

fs is strictly positive in (— cxd, a„ ((5)) U (a((5), a+((5)), (37) 
fs is strictly negative in (a ^ (5), a{6)) L) {a ^{5), +00). 

• Set 

M5) := f^{a{6)) = na{6)), 

then there exists a positive constant, which we denote again by C, such 
that 

\^i{^)-^^\<C\S\. (3.8) 

In order to construct a pair of sub and super-solutions for Problem (P^) 
we define Y{t, ^; 5) as the solution of the ordinary differential equation 

( Yr{T,^;6) =fs{Y{T,^;6)) for r > 0, 

\ y(o,e;<5) =e, 

for 5 £ {—5o,5o) and ^ G (— 2Co,2Co), where Cq has been chosen in (II. 2p . 
We present below basic properties of Y. 
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Lemma 3.3. We have > 0, for all ^ G {-2Co,2Co)\{a-{6),a{6),a+{6)}, 
all 6 G (— (5o,(5o) and all r > 0. Furthermore, 



f5iY{T,^;5)) 

Proof. We differentiate (13. 9p witli respect to ^ to obtain 

r = Yi:f'{Y), 

\ Y^i0,^;5) = l, 
wliich is integrated as follows: 



y5(r,C;<5) =exp / f'{Y{s,C;S))ds 

L Jo 



> 0. 

Then differentiating (13. 9p with respect to r, we obtain 

f Yrr = Yrr{Y), 

\ Yr{0,^;6)=fs{O, 

which in turn implies 

Yr{T,C;5)=fs{0^W \ r f'{Y{s,i;5))ds 

L JO 

which enables to conclude. 

We define a function A(r, ^; 5) by 



(3.10) 



□ 



nY{T,^-5))-no 

fsiO 



(3.11) 



Lemma 3.4. We have, for all £, £ (-2Co,2Co) \ {a^{S),a{6),a+{6)}, all 
6 £ {—So, So) and all r > 0, 



A{t,^;6)= / f"{Y{s,^;6))Y^{s,^;6)ds. 
Jo 

Proof. We differentiate the equality of Lemma 13.31 with respect to ^ to 
obtain 

Y^^{t,C;S)=A{t,C;S)Y^{t,^;6). (3.12) 
Then differentiating (13.101) with respect to ^ yields 



□ 



n? = ^« / f"{Y{s,^;6))Y^{s,^;6)ds. 
Jo 

These two last results complete the proof of Lemma 13.41 



Next we prove estimates on the growth of Y, A and theirs derivatives. 
We first consider the case where the initial value ^ is far from the stable 
equilibria, more precisely when it lies between rj and 1 — rj. 
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Lemma 3.5. Let r] be arbitrary. Then there exist positive constants 6q = 
6o{r]), Ci = Ci{r]), C2 = C2iv) "^^^c? C3 = 6*3(7/) such that, for all 6 G 
{-60,60), for all r > 0, 

• if i ^ {o,{6), 1 — rf) then, for every r > such that Y[t,^;5) remains 
in the interval {a{5), 1 — rj), we have 

C'^eM('5)r < Y/_{t, 6) < Cie^^^^^^ (3.13) 

and 

|A(r,e;<5)| <C3(e^(^)--l); (3.14) 

• if S, ^ iv^'^i^)) then, for every r > such that Y{t,S,;S) remains in 
the interval {rj,a{6)), ()3.13p and (I3.14|) hold as well. 

Proof. We take ^ G {a{6), 1 — ?/) and suppose that for s G (0, r), Y{s,S,; S) 
remains in the interval {a{6), 1 — rj). Integrating the equahty 

Yr{s,C;6) _^ 
fs{Yis,^;6)) 

from to r and using the change of variable q = Y{s,(,; 6) leads to 
Moreover, the equality in Lemma 13.31 enables to write 

r.y(r,C;<5) rn 



fsiq) fs{q) 



{r,t,;0) 

h{6)t + j hs{q)dq, 



where 

fsiq) 

In view of (|3.8p . respectively (j3.6|) . we can choose 60 = 6o{i]) > small 
enough so that, for all 6 G [— (^O)'^o]) we have fi{6) > n/2 > 0, respectively 
{a{6),l — rj] C {a{6), a+{6)). Since 

f's'iai6)) r{a{5)) 

we see that the function (q, 6) ^ hs{q) is continuous in the compact region 
{ l^^l ^ 60, a{5) < q < 1 — r]}. It follows that |/i5(Q')| is bounded by a constant 
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H = H{r]) as {q, 5) varies in this region. Since |y(r, 6) — ^| takes its values 
in the interval [0,1 — rj — a{5)] C [0, 1], it follows from (j3.16p that 

/u((5)r -H< InY^iT, 5) < fi{6)T + H, 

which, in turn, proves (|3.13|) . Next Lemma 13.41 and (|3.13|) yield 

\A{t,^;6)\ <\\n\L^io,i) rC2e^'^'^'ds 

Jo 

- KS) ^ ^ 
<-|irilL-(o,i)C'2(e^(')^-l), 

which completes the proof of ()3.14p . The case where ^ and Y{t, ^; 5) are in 
{r],a{S)) is similar and omitted. □ 

Corollary 3.6. Let r] be arbitrary. Then there exist positive constants 6q = 
6o{r]), Ci = Ciirf) and C2 = C2{r]) such that, for all 5 S {—5o,5o), for all 
T > 0, 

• if £ {a{S), 1 — rj) then, for every r > such that Y{t,S^;6) remains 
in the interval {a{S), 1 — rj), we have 

C7^eM5)r(^ - a{6)) < Y{t, S) - a{6) < C2e'^(^)"(^ - 0(6)), (3.17) 



if ^ £ {ri,a{6)) then, for every r > such that Y{t,£,;6) remains in 
the interval {rj,a{6)), we have 

C2e''^^'>^{^ - a{S)) < Y{t, S) - a{S) < de'^'^^^^i^ - a{6)). (3.18) 



Proof. In view of ()3.8p . respectively ()3.6|) . we can choose 60 = 6o{r]) > 
small enough so that, for all 6 G [— (5o,(5o], we have fi{6) > n/2 > 0, 
respectively {a{6),l — 1]] C {a{6),a+{6)). Since 

■^^^'^^ ^1(6) as g— >a(5), 



q - a{6) 



it follows that {q, 5) ^ f5{(l)/{<l — «(^)) is a strictly positive and continuous 
function in the compact region {\5\ <5q, a{5) < q < 1 — ?/}, which insures 
the existence of constants Bi = Bi[rf) > and B2 = B2{r}) > such that, 
for all q G (a((5), 1 — ??), all 5 G {—5q, 5q), 

Bi{q - a{6)) < fs{q) < B2{q - a{5)). (3.19) 
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We write the inequalities ()3.19p for q = Y{t,S,;6) G 1 — f]) and then 

for g = ^ G ('^C'^)) 1 ~ '^)) which, together with Lemma |3.3^ imphes that 

^{Y{T,^;S)-a{6)) < {^-a{S))Y^{T,^;S) < ^{Y{r,^;S) - a{6)). 

-D2 Bi 

In view of (jS.lSp . this completes the proof of inequalities ()3.17p . The proof 
of (|3.18p is similar and omitted. □ 

We now present estimates in the case that the initial value ^ is smaller 
than T] or larger than 1 — i]. 

Lemma 3.7. Let rj and M > be arbitrary. Then there exist positive 
constants 5o = 5o{t],M) and C4 = Ci{M) such that, for all 5 G {—5o,6o), 

• if C G [1 — ?7, 1 + M], then, for all t > 0, Y{t,£,]6) remains in the 
interval [1 — 77, 1 + M] and 

\A{t,^;6)\ <CiT for T>0; (3.20) 

* if G [—M,r]], then, for all r > 0, Y{t,£,;5) remains in the interval 
[—M, rj] and (|3.2Up holds as well. 

Proof. Since the two statements can be treated in the same way, we will 
only prove the former. The fact that Y{t, ^; 5), the solution of the ordinary 
differential equation (j3.9p . remains in the interval [1 — r/, 1 + M] directly 
follows from the bistable properties of fs, or, more precisely, from the sign 
conditions fs{l - r]) > 0, fs{l + M) < valid if 60 = 6o{r],M) is smaU 
enough. 

To prove (|3.20p . suppose first that ^ G [a+{6),l + M]. By the above 
arguments, Y(t, 6) remains in this interval. Moreover /' is negative in 
this interval. Hence, it follows from (|3.10p that Y^{t,(,;6) < 1. We then use 
Lemma 13.41 to deduce that 

\A{t,^;6)\ < \\f"\\L^(-M,i+M)T =: C4T. 

The case ^ G [1 — r/, a_|_((5)] being similar, this completes the proof of the 
lemma. □ 

Now we choose the constant M in the above lemma sufficiently large 
so that [-2Co,2Co] C [-M, 1 + M], and fix M hereafter. Therefore the 
constant C4 is fixed as well. Using the fact that r 1-^ T{e^^^^'^ — is 
uniformly bounded for S G {—6o,So), with 60 small enough (see ()3.8p ). and 
for r > 0, one can easily deduce from (|3.14p and (|3.2Up the following general 
estimate. 

Lemma 3.8. Let rj be arbitrary and let Co be the constant defined in ()1.2p . 
Then there exist positive constants 60 = So{r]), C5 = 6*5(77) such that, for all 
5 G {-So, 60), all C G (-2Co, 2Co) and all t > 0, 

|A(r,e;<5)| <C5(e^('^)--l). 
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3.2 Construction of sub and super-solutions 

We now use Y to construct a pair of sub- and super-solutions for the proof 
of the generation of interface theorem. We set 



(x, t) = y ( — , uo(x) ± eV(±eG, ±eG] , 



We 

where the constant G is defined by 

G= sup \g{u)\, 

«e[-2Co,2Co] 

and the function r((5, r) is given by 

r(5,r)=C6(e^(^)^-l). 
For simpHcity, we make the fohowing additional assumption: 

dv 



(3.21) 



on do.. 



(3.22) 



In the general case where ()3.22p does not necessary hold, we have to slightly 
modify near the boundary 50. This will be discussed in the next remark. 

Lemma 3.9. There exist positive constants Eq and Cg such that for all e £ 
(0,eo), the functions w~ and wf are respectively sub- and super-solutions 
for Problem (-P^), in the domain 

{{x,t) gQt, x£n, 0<t< lne|}. 

Proof. First, (|3.22|) implies the homogeneous Neumann boundary condition 

^ - on d^lx (0,+oo). 



du 



Let v'^ be such that 



V Of 

According to Definition 12.21 what we have to show is 



(3.23) 



LAwt] ■■= {wt)t - Awt + V • {wtVxiv')) - ^fe{wt) > 0. 

Let Cg be a positive constant which does not depend on e. If eq is suffi- 
ciently small, we note that iteG G (— 5oi ^o) ^-nd that, in the range < i < 



|e2C7g(e/^(±^«)*/=' - 1)1 < e2Q(g-/.(±eG)//. _ 1) < 
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which imphes that 



uo{x)±e\{±eG,t/e^) G (-2Co,2Co)- 

These observations allow us to use the results of the previous subsection 
with r = t/e'^, ^ = uo{x) + e^r(eG, t/e^) and 6 = eG. In particular, setting 
Fi := ||/'||l°°(-2Co,2Co)> this implies, using (|3.10!) . that 



< < e 



Straightforward computations yield 

and then, in view of the ordinary differential equation (|3.9|) . playing the 
role of (5, 



+ V^zo • Vx(^;^) + ^Ax(t;")J . (3.24) 

By the definition of G the first term above is positive. Now, using the choice 
of Go in (jl.2p . the fact that Y^^jY^ = ^ and Lemma [3^ we obtain, for a 
G5 independent of e, 

[«^^ ] > y^ [G6Ai(eG)e'^(^G)*/"' - Go - CgCe'^^"^)*/^' - IjGo^ 
-Go|Vx(T;^)|-2Goe^i^|Ax(t'^)| . 

Moreover, the inequalities in (13.230 can be written as — A?;^ + 7f^ = /i^, with 
— 2Go <hf< 2Go, so that the standard theory of elliptic equations gives a 
uniform bound M for |Vf^| and |Au^|. Hence, using the smoothness of 
X, we have a uniform bound M' for |Vx(f'^)| and |Ax(w^)|. It follows that 

L,^ [wt] > n [(G6Ai(eG) - G5Go2)e^(^^)*/^' - Go+GgGo' - Gq Af - 2Goe^^'^M' 

Hence, in view of (j3.8p . we have, for eq small enough (recall that 1^ > 0), 



Lwt > Y^ 



{Ce^l^ - C^Ci) - Go - GqM' - 2Goe^i^ M' 



>0, 



for Cq large enough, so that wf is a super-solution for Problem (P^)- We 
omit the proof that wZ^ is a sub-solution. □ 
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Now, since wf{x,0) = Y{0,uo{x);±eG) = uq{x), the comparison prin- 
ciple set in Proposition 12.31 asserts that, for ah x G il, for aU < t < 

w~{x,t) <u''{x,t) <wf{x,t). (3.25) 

Remark 3.10. In the more general case where ()3.22p is not valid, one can 
proceed in the following way: in view of ()1.3p and (II. 4p there exist positive 
constants di and p such that uo{x) < 1/2 — p if d{x,dfl) < di. Let x 
be a smooth cut-off function defined on [0,-|-oo) such that < x ^ 1) 
x(0) = x'(0) = and x(z) = 1 for z > di. Then define 

u+{x) = x(.d{x,dn))uo{x) + (1 - x(.d{x,dn)){l/2 - p), 
Uq (x) = xi^ix, (?r2))no(x) + (1 — xid-ix, dO,)) minMo(a;)- 
Clearly, Uq < uq < Uq , and both Uq satisfy (|3.22p . Now we set 
wf{x,t)=Y(^^, u^{x) ± e\{±eG, ^y,±eGy 

Then the same argument as in Lemma 13.91 shows that {w~,wf) is a pair of 
sub and super-solutions for Problem (P^)- Furthermore, since 'w~(x,0) = 
Uq{x) < uo{x) < Uq{x) = 'wf{x,0), Proposition 12.31 asserts that, for all 
X G SI, for all < t < /U~^e^| lne|, we have w^{x, t) < u^{x, t) < wf{x, t). □ 

3.3 Proof of Theorem 

In order to prove Theorem l3.1l we first present a key estimate on the function 
Y after a time of order r ~ | lne|. 

Lemma 3.11. Let r] be arbitrary; there exist positive constants Eq = £o{r]) 
and Cy = Ci^rj) such that, for all e G (0, eo)? 

. for aine{-2Co,2Co), 

-r] <Y{p,^^\lne\,(,;±eG) <l + r], (3.26) 

• for all G (— 2Co,2Co) such that |^ — ^| > Cje, we have that 

if ^>l/2 + C7e then Y{p'^\\ne\,^;±eG) > 1 - rj, (3.27) 
if i<\l2-C-je then Y {pr^\\^e\,i\±eG) <r]. (3.28) 

Proof. We first prove (|3.27p . In view of (j3.6p . we have, for C7 large enough, 
l/2-|-C7e > a(±eG) -I- ^Cye, for all e G (0, eo); with eq smaU enough. Hence 
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for ^ > 1/2 + 6*76, as long as y(r, ^; iteG) has not reached l — rj, we can use 
p.l7p to deduce that 

Y{t, ±eG) > a{±eG) + Cie^(±^^)^(^ - a{±eG)) 

> a{±eG) + ICiGjee''^^^^^^ 

> I - eGG + iCiCyee'^^^^^)^ 

> 1 - 77 

provided that 

, 1 , 1/2 - 7? + CGe 
T > T := In . 

fi{±eG) GiG7e/2 

To complete the proof of (j3.27p we must choose C7 so that lne| — > 0. 
A simple computation shows that 

Thanks to (|3.8|) . as e — > 0, the first term above is of order e\ lne| and the 
second one of order 1. Hence, for C7 large enough, the quantity lne| — 
is positive, for all e G {0,eo), with £q small enough. The proof of (|3.28|) is 
similar and omitted. 

Next we prove (|3.26p . First note that, by taking eo small enough, the 
stable equilibria of f±eG^ namely Q;_(±eG) and a+(±eG), are in [— r/, l + rj]. 
Hence, f±eG being a bistable function, if we leave from a ^ G [—rj, 1 + r]] 
then y(r, ibeG) will remain in the interval [— 1 + f?]. Now suppose that 
1 + 77 < ^ < 2Co (note that this work is useless if 2Go < 1 + ??) • We check 
below that ¥{^-^1 lne|, ^; ±eG) <l + r). As long asl + r]<Y < 2Go, <^ 
leads to the inequality Y-j- < /(I + t]) + eG < ^/(l + ??) < 0, for eq = eo(?7) 
small enough. By integration from to r, it follows that 

YiT,C;±eG) <e + i/(l + 7?)r 

< 2Go + i/(l + 7?)r 

< 1 + r?, 

provided that 

2Co - 1 - r? 

^- -/(l + 7?)/2' 

and a fortiori for r = lne|, which completes the proof of (I3.26h . □ 
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We are now ready to prove Theorem 13.11 By setting t = fi lne| in 
(|3.25p . we obtain 



Yyfi ^\lne\, uo{x) — e'^r{—eG, iJ, ^\lne\); —eG ^ 

< M^(x,^i"^e2|lne|) < Y (^fx~'^\lne\,uo{x) + e^r(eG, ^"^| In e|); . 

(3.29) 



In view of (IHI 



lim^i^i^i±^lne = 0, (3.30) 
SO that, for eq small enough, we have 

eV(±eG,^-i|lne|) =C6e(e(^-^(±^^«/^ -e) G {^Gee,^Gee). 

It follows that uo{x) ± £'^r{±£G, fi~^\lne\) G (— 2Co,2Co). Hence the result 
()3.2p of Theorem 13.11 is a direct consequence of ()3.26p and (13.290 . 

Next we prove (j3.3p . We take x G such that uo{x) > 1/2 + Mqe so 
that 

uoix) - e'^r{-eG,fi~^\lne\) > 1/2 + Mqe - ICqE 

> 1/2 + Cre, 

if we choose Mq large enough. Using (j3.29p and (j3.27p we obtain (|3.3p . which 
completes the proof of Theorem 13.11 □ 



4 Motion of interface 



We have seen in Section [3] that, after a very short time, the solution 
develops a clear transition layer. In the present section, we show that it 
persists and that its law of motion is well approximated by the interface 
equation in (P") obtained by formal asymptotic expansions in subsection 

EH 

More precisely, take the first term of the formal asymptotic expansion 
(1231) formal expansion of the solution: 

u%x,t) « u-{x,t):=Uo[^^), (4.1) 

where Uq is defined in (|2.6p . The right-hand side is a function having a 
well-developed transition layer, and its interface lies exactly on Ft- We show 
that this function is a good approximation of the solution; more precisely: 

If u'^ becomes very close to at some time moment t = to, then 
it stays close to v!^ for the rest of time. Consequently, Ff evolves 
roughly like Tt- 
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To that purpose, we will construct a pair of sub- and super-solutions 
and uf for Problem (P^) by slightly modifying u^. It then follows that, if 
the solution satisfies 

u~{x,to) < u''{x,to) < uf{x,to), 

for some Iq >0, then 

uj{x,t) < u^{x,t) < uf{x,t), 

for to < t < T. As a result, since both uf,u~ stay close to u^, the solution 
also stays close to for to <t <T. 

4.1 Construction of sub- and super-solutions 

To begin with we present mathematical tools which are essential for the 
construction of sub and super-solutions. 

A modified signed distance function. Rather than working with the 
usual signed distance function d, defined in (12. ip . we define a "cut-off signed 
distance function" d as follows. Choose do > small enough so that d{-, •) 
is smooth in the tubular neighborhood of F 

G Q^, \d{x,t)\ <3do}, 

and such that 

dist(rt, dn) > 4do for all t G [0, T]. (4.2) 
Next let C{s) be a smooth increasing function on M such that 

s if \s\ < 2do 

C(s) = { — 3(io if s < —3do 
3do if s > 3do- 

We define the cut-off signed distance function d by 

d{x,t) = C{d{x,t)). (4.3) 

Note that \Vd\ = 1 in the region {{x,t) G Qt, \d{x,t)\ < 2do} and that, in 
view of the above definition, Vd = in a neighborhood of dil.. Note also 
that the equation of motion interface in (P^), which is equivalent to ()2.9p . 
is now written as 

dt = Ad-Vd- Vx(f °) - V2a on Ft. (4.4) 

Construction. We look for a pair of sub- and super-solutions for Prob- 
lem (P^) of the form 

„f(,,,) = [;„(*!i^)p*))±,(,), (4.5) 
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where Uq is the solution of (12. 6p . and where 

p{t) = -e-^Ve' + e-^* + K, 

(4.6) 

Note that q = ae^pt- Let us remark that the construction (j4.5p is more 
precise than the procedure of only taking a zeroth order term of the form 
C/q, since we have shown in the formal derivation that the first order term 
C/i in ()2.3p vanishes. It is clear from the definition of uf that 

, f 1 for all (x, t) G Q^}) 

]iuiut{x,t) = { 7 (4.7) 

i for all {x,t) G Q^"\ 

The main result of this section is the following. 

Lemma 4.1. There exist positive constants /3, a with the following proper- 
ties. For any K > I, we can find positive constants Eq and L such that, for 
any e G (0, eo); the functions u~ and uf are respectively sub- and super- 
solutions for Problem (P^) in the range x£O,,0<t<T. 

4.2 Proof of Lemma 14.11 

First, since Vd = in a neighborhood of dil., we have the homogeneous 
Neumann boundary condition 

= on dQx [0,T]. 

Let v'^ be such that (12. lip holds. We have to show that 

L,4ut] := iut)t - Aut + Vut ■ Vx{v') + ufAxiv') - ^/.(n+) > 0, 

the proof of inequality Ljjs[u^] < following by the same arguments. 

Computation of LyE[uf]. By straightforward computations we obtain the 
following terms: 

iut)t = Uo'{^-pt) + qt, 

Vut = Uo'^, 

AuJ = Uq 5 h Uq , 

where the function Uq, as well as its derivatives, is taken at the point 
[d(x,t) — ep{t))/e. We also use expansions of the reaction terms: 

f{ut) = f{UQ) + qf'{UQ) + ^q'f"{e), 

g{ut)=g{UQ)+qg'{uj), 
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where 6{x,t) and uj{x,t) are some functions satisfying Uq < 9 < uf , Uq < 
uj < uf. Combining the above expressions with equation (j2.6p and the fact 
that V2aUo' + g{UQ) = 0, we obtain 

L^.[n+] =Ei + --- + E5, 

where: 



El = -^q[f'{Uo) + Iqf'ie)] - Uo'pt + qt, 



E, = %{l-\Vd\% 

E3 = ^{dt- Ad + Vd- Vxiv°) + V2a), 

E4 = -^qg'iuj), 

E, = ^Vd ■ Vixiv') - x(f°)) + utAxiv')- 

In order to estimate the terms above, we first present some useful in- 
equahties. As /'(O) = /'(I) = —1/2, we can find strictly positive constants 
h and m such that 

if Uo{z) G [0, 6] U [1 - 6, 1] then f {Uq{z)) < -m. (4.8) 

On the other hand, since the region {z E M | Uq{z) G [6, 1 — 6] } is compact 
and since Uq < on M, there exists a constant oi > such that 

if Uq{z) G [6, 1 - h] then Uq{z) < -ai. (4.9) 

We then define 

F= sup \fiz)\ + \f'iz)\ + \f"{z)\, (4.10) 

-l<z<2 

/3 = J, (4.11) 

and choose a that satisfies 

< cr < min (do, 0-1,(12), (4.12) 

where 

ai 1 4/3 

•^0 := —p^, := -7——, 0-2 : = 



m + F' /3+1' F(/? + l)' 

Hence, combining (|4.8|) and (|4.9|) . we obtain, using that a < ao, 

-Uo'iz)-af'{Uo{z))>4a/3 for - 00 < z < 00. (4.13) 
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Now let K > 1 be arbitrary. In what follows we will show that L^e [uf] > 
provided that the constants Eq and L are appropriately chosen. From now 
on, we suppose that the following inequality is satisfied: 



egLe^' < 1 . (4.14) 

Then, given any e G (0, Eq)) since a < ai, we have < q{t) < 1, hence, 
recahing that < Uq < 1, 

-1 <uf{x,t) <2. (4.15) 



We first estimate the term Ei. A direct computation gives 
El = ^ e-f^'''\l - al3) + Le^\l + E^aL), 



I = -Uo' - afm - ^f"{e){(3e-^'/^' + e^Le^*). 



where 

2 

In virtue of (j4.13p and (|4.15|) . we obtain 

/ > 4a/3 - Yi^(/? + e'Le^ 
Then, in view of (|4.14|) . using that cj < cj2, we have 

/ > 2a(3. 

Consequently, the following inequality holds. 

The term £'2- First, in the points where where |(i(a;,t)| < do, we have that 
\Vd\ = 1 so that E2 = 0. Next we consider the points where |d(x,t)| > do. 
We deduce from Lemma l2.1l that: 

\E2\ < -5(l + ||Vd||^)e-^l'^+^^'l/^ 
C 

In view of the definition of p in (j4.6p . we have that 0<A'— 1 < p < e^^ + K, 
and suppose from now that the following assumption holds: 

e^^ + K<A. (4.16) 

2eo 



<^(l + l|Vd||L)e-"("°/^-l^l). 
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Then — — \p\ > —, so that 

£ 2e 
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<C2:=-^(l + ||Vd||L). 
(eAdo)"' 



Next we consider the term £^3. We recall that 

(ij - Ad + Vd • Vx(i'°) + \/2a = on Ft = {x G J^, d(x, t) = 0}. 

Since is of class C^^^'—, for any t?' G (0,1), and since the interface 
Ft is of class C^^'''^~, the functions Vd, Ad, dt and Vx(t''') are Lipschitz 
continuous near Vt- It then follows, from the mean value theorem applied 
separately on both sides of F^, that there exists Nq> ^ such that: 

|(dj- Ad + Vd-Vx(i'°) + \/2a)(x,t)| <iVo|d(x,t)| for all (x,t)GQr. 

Applying Lemma 12.11 we deduce that 

|^3| < jVgC'M^il)le-A|rf(x,t)/e+p(j)| 
e 

< AfoCmaXj^eM \y\er^\yMAt)\ 

< AoCmax(|p(t)|,i) 

Taking the expression of p into account, we see that < e^* + A', which 

implies 

lii^sl <C3(e^* + /^) + C3', 
where C3 := AqC and C3' := AqC/A. 

The term £^4. We set G\ := \\g'\\L°°{-i,2) and, substituting the expression 
for q, obtain that 

\E,\ <aGi(^e-^*/^'+eLe^*) 

e 

We continue with the term £5. This term requires a more delicate 
analysis. We need a precise estimate of t;^ — . We recall that satisfies 
—A?;'' + = , with a step function discontinuous when crossing the 
interface. 
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Lemma 4.2. There exists a positive constant Cq such that, for all (x, t) G 
Qt, 

{\v^\ + \Vv^\ + |Av=|)(x,t) < Cg, (4.17) 
{\v' - + |Vd • y{v' - {x,t) < Cciepit) + q{t)). (4.18) 

We postpone the proof of this lemma and pursue the proof of Lemma 14.11 
Using the smoothness of x ()4.17p . we obtain a uniform bound Cg' for 
Ax{v^)- Moreover, we write 

Vd • V{x{v') - xiv'^)) = X{v')^d ■ V{v' - v'>) + {x'{v') - x'{v'))^d ■ Vv'. 

(4.19) 

Since is of class C^~^^ , for any -d' S (0, 1), there exists a constant, 
which we denote again by Cq, such that 

II^°IIl°°(Qt) + II^^°IIl°°(Qt) ^ ^G, 
which, combined with (j4.19p . yields 

\Vd-V{x{v')-xiv'))\ < \\x'\\oo\Vd-Viv'-v^)\ + \v'-v'\ ||x"l|oo || VdHooCc, 

where the L°°-norms of x' c^nd x" are considered on the interval {—Cg, Cg)- 
It follows from the above inequality and ()4.18p that there exists a constant 
Cg" such that, for all {x,t) G Qt, 

\Vd.V{xiv^) - xiv'))\{x,t) < CG"{ep{t)+q{t)). 

Hence, using ()4.15p and the above estimates, we obtain, 

|^5| < ^CG"{ep{t)+q{t)) + 2CG'. 

Then, substituting the expressions for p and g, we easily obtain positive 
constants C5, C5' and C5" such that 

m < C5 + ^e-W^' + C5"(l + eL)e^*. 

Completion of the proof. Collecting the above estimates of Ei — 
yields 

LAut] > ~ " '^'' e-f''/^' + {L{C,'-eC,'-eC,")-Cs-C,")e^'-C,, 
where C7 := C2 + KC3 + C3 + C5. Now, we set 

T 4eo' 
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which, for eq small enough, validates assumptions ()4.14p and (I4.16p . If Eq 
is chosen sufficiently small (i.e. L large enough), Ci/e^ — (C4 + Cs')/^ is 
positive, Ci - eC^ - eCg" > \Ci , and 

L.^ [ut] > i^LCi' -C3- C5"] e^* - Cr 

> \LCi' - C7 

> 0. 

The proof of Lemma |4. H is now completed, with the choice of the constants 
/3,cr as in gH]), ([112]). □ 



4.3 Proof of Lemma 14.21 

Lemma 14.21 is inspired by Lemma 4.9 in [1]. Since our pair of sub- and 
super-solutions is different from the one in we need to perform some 
minor changes. First we give a useful estimate on "shifted Uq\ 

Lemma 4.3. For all a G M, all z ^ we have 

\Uo{z + a) - X]-oo,o](^)l < Ce-^l^+'^l + X[-a,a](^) 

Proof. Let us give the proof for a > 0. We distinguish three cases and use 
the estimates of Lemma l2.1[ For z < —a, we have |C/o(-2+fl)~l| ^ Ce""^!^"*""!. 
For -a < z < 0, we have \Uo{z + a) - l\ < \Uo{z + a)| + 1 < Ce-^l^+^l + 1. 
For z > 0, we have |C/o(-2^ + ^ Ce~'^l^"'"'*L We proceed in the same way 
for a < 0. □ 

We turn to the proof of Lemma 14. 2i First, we recall that is such 
that ()2.1ip holds; hence, in view of ()4.15p . the estimate (I4.17P is a direct 
consequence of the standard theory of elliptic equations. Next we prove 
(j4.18p . The function w = ■= v'^ — is solution of 

—Aw + '^w = h on Qt, 

dw (^-20) 
— = on5Ox(0,r), 
ov 

with u~ — < h = h'^ < uf — u'^, where is the step function defined 
by u^{x,t) = X{d{x,t)<o}- The key idea of the proof is the fact that h is 
exponentially small with respect to e, except possibly in a thin neighborhood 
of Tt of width of order ep{t). More precisely, from the definitions of in 
()4.5p and from the above lemma for z = d{x, t)/e and a = zizp{t), we deduce 
that 

\hix,t)\ < C(e-^l'^(-'*)/^+PWI +e-^l'^(-'*)/^-*'WI) +X{M(.,t)|<.pW} 

(4.21) 
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By linearity, we successively consider equation ()4.20p with the various terms 
appearing in the right-hand side of (j4.2ip . By the standard elliptic esti- 
mates, the solution w of ()4.20p satisfies 

\w{x,t)\ +\Vw{x,t)\ < C sup \h{y,t)\, (4.22) 

yen 

which gives the term Ccqit) that appears in the right-hand side of inequality 
(|4.18p for h{y,t) = q{t). We now suppose that the function h satisfies one 
of the three following assumptions: 

(H^) \h{y,t)\<X{\d^y^t)\<ep{t)} 
{Hf) <exp(-A|^^±p(i)|), 



and write 



h{y,t) = h{y,t)x{\diy,t)\<do} + Ky,'t)X{\d[y,t)\>do}- 



We first consider the term h{y,t)x{\d{y,t)\>do}- virtue of (j4.16p . we 
have 

< K -1 <p{t) <do/2eo. (4.23) 

Under assumption (Hi), it follows that h is supported in {\d{y,t)\ < do/2}, 
which implies h{y,t)x{\d{y,t)\>do} = 0- Moreover, under assumption (H^), 
using again (j4.23p . 

\h{y, t)\X{\d{y,t)\>do} < exp [ - X{do/e - p{t))] 
< exp(— A(io/2e) 

^ 2 
~ Xdoe 

Thus, under either of the assumptions (Hi) or (H^), the estimate ()4.18p — 
for the term h{y,t)x{\d(y,t)\>do} — directly follows from inequality (|4.22p . 

From now on, we assume that h is supported in {|(i(?/, t)| < do}. We 
have that 



and 



w{x,t)= G{x,y)h{y,t)dy, 

J\d(y,t)\<do 

Vd{x,t)-Vw{x,t)= / {V^G{x,y) ■Vdix,t))h{y,t)dy, 

J\d{y,t)\<do 

where G is the Green's function associated to the homogeneous Neumann 
boundary value problem on Vl for the operator — A -|- 7. More precisely. 
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G{x,y) = g-y{\x — y\) + H-y{x,y), where g'y{\x — y\) is the Green's function 
associated to the operator —A + 7 on and where Hy{x, y) is smooth for 
X and y far away from dVl. It is known that g^/ is the Bessel function defined 

by 

g^[r) = CN e 2se ^2s 2 

./n s 



with Cat > a normahzation constant. We use the following estimates (see 
01): 

for N >3 



\G{x,y)\ < { 



y — x\^ ^ 
_C| ln|y - x\ 



|V..G(.,,).Vd(x,t)|< ^'"(^'^)-5"'^)' + 



for iV = 2, 
C 



\y -X 



N 



\y-x 



N-2 



(4.24) 

for N>2. 

(4.25) 



This last inequality follows from 



\V^G{x,y)-Vd{x,t)\ < 



C\Vdix,t) ■ {y - x)\ 



\y -X 



N 



and from d{y,t) — d{x,t) = Vd{x,t) ■ {y — x) + 0{\y — xp). Now, under 
respectively assumptions {Hi), (H^), we define a function h = on R x 
[0,T], respectively by 



h{r,t) := < 



H\r\ < epit)} 
exp( - X\-±p{t)\ 



(4.26) 



Note that \h{y,t)\ < h{d{y,t),t). Moreover, using ()4.23p . straightforward 
computations show that, under either of the assumptions (Hi) or (H^), 
there exists C > such that 



< 



d-O 



do 



h{r,t)dr < Cep{t), 



(4.27) 



which is an analogue of (4.20) in [1] . The end of the proof is now identical to 
that of Lemma 4.10 in [3]. We omit the details and refer to this article. □ 



5 Proof of the main results 

In this section, we prove our main results by fitting the two pairs of sub- and 
super-solutions, constructed for the study of the generation and the motion 
of interface, into each other. 
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5.1 Proof of Theorem \TM 

Let 1] G (0, 1/4) be arbitrary. Choose (3 and a that satisfy ()4.1ip . (14.121) and 

a/3<|. (5.1) 

By the generation of interface Theorem 13. H there exist positive constants 
eo and Mq such that ()3.2p . ()3.3p and (13. 4p hold with the constant r] replaced 
by a (3/2. Since Vuq • n / everywhere on the initial interface Tq = {x £ 
Q., uo{x) = 1/2} and since Tq is a compact hypersurface, we can find a 
positive constant Mi such that 

if do{x)> Mie then no(3;) < 1/2 - Mq^, 
if do{x)<-Mie then no(2;) > 1/2 + Moe. 

Here do{x) := (i(x,0) denotes the cut-off signed distance function associated 
with the hypersurface Tq. Now we define functions {x), H~ (x) by 



H+{x) 



H-{x) 



1 + if do{x) < Mie 

\a(3 if dQ{x) > Mie, 

1 - icj/3 if do(x) < -Mie 

-\(t(3 if do(x) > -Mie. 



Then from the above observation we see that 

H^{x) < u^ix, fi-^e"^ \lne\) < H+{x) for x e n. (5.3) 
Next we fix a sufficiently large constant K > 1 such that 

Uo{Mi -K)>1-^ and Uo{-Mi + K) < ^. (5.4) 

For this K, we choose eg and L as in Lemma l4.ll We claim that 

u-{x,0) < H-{x), if+(x) < ti+(x,0) for x £ n. (5.5) 

We only prove the former inequality, as the proof of the latter is virtually 
the same. Then it amounts to showing that 

uj{x,0) = Uo{^^ + K) -a{p + e^L) < H~{x). (5.6) 

In the range where do{x) > —Mie, the second inequality in (15. 4p and the 
fact that C/q is a decreasing function imply 

Uo{^^ + K)- a{p + e^L) < \ap - aP < H-{x). 
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On the other hand, in the range where do(x) < —Mie, we have 
Uo{^^ + K)-a{(3 + e^L)< l-a(3 < H-{x). 



This proves ()5.6p . hence (15. 5p is estabhshed. 
Combining (j5.3p and (jS.Sp . we obtain 

(a;,0) < n'^(x,;U~-^e^| lne|) < u+(x,0). 

Since u~ and are sub- and super-solutions for Problem (/"^) thanks to 
Lemma l4.lt the comparison principle yields 

uj{x,t) <u^x,t + f) <uf{x,t) for 0<^<^-^^ (5.7) 

where = fi~^e'^\lne\. Note that, in view of (|4.7p . this is enough to prove 
Corollary 11.41 Now let C be a positive constant such that 

C/o(-C + e^^ + i^) > 1- I and Uo{C - e^^ - K) < ^. (5.8) 

One then easily checks, using successively (15. 7p . (14. 5p . ()5.8p and (15. ip . that, 
for eo small enough, for < t < T — t^, we have 

if d{x,t)> Ce then ^"^(2;, t -h i^) < ?/, 
if d{x,t)<-C£ then -(/^(x, t + i^) > 1 - r/, 

and 

ti''(x,t-Ft^) e [-r/, 1 -Fry], 
which completes the proof of Theorem II. 3[ □ 



5.2 Proof of Theorem [T3] 

In the case where /i~^e^|lne| < t < T, the assertion of the theorem is a 
direct consequence of Theorem 1 1.31 All we have to consider is the case where 
0<t< n'^e^l lne\. We shah use the sub- and super-solutions constructed 
for the study of the generation of interface in Section EJ To that purpose, 
we first prove the following lemma concerning Y{t, ^; 6), the solution of the 
ordinary differential equation (j3.9p . in the initial time interval. 

Lemma 5.1. There exist constants Cg > and £q > such that, for all 
e e (0,eo), 

if (,>l/2 + Cs£ then y(r, ^; ±eG) > 1/2 for < r < ^ti^^l lne|, 

if C<l/2-C8e then Y{t,C;±£G) < 1/2 for < t < fi'^\lne\. 

(5.10) 
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Proof. We only prove the first inequality. Assume ^ > l/2 + C8e. By (13. 6p . 
for Cs > CG, we have that ^ > 1/2 + Cse > a{±eG). It then follows from 
([3X7]) that 

Y{t, e; ±eG) > a{±eG) + Cie'^(±=^)^(l/2 + Cge - a{±eG)) 

> 1/2 - CGe + Ci(-CGe + Cge) 

> l/2 + e(CiC8-CG(Ci + l)) 

> 1/2, 

provided that Cs is sufficiently large. □ 

Now we turn to the proof of Theorem II. 5i We first claim that there 
exists a positive constant M2 such that for all t G [0,fi~^e'^ \ lne|], 

TtcMM2e(Xo). (5.11) 

To see this, we choose Mq' large enough, so that Mq > Cs + 2Cq, where 
Cq is as in Lemma 13.91 As is done for (|5.2|) . there is a positive constant M2 
such that 

if do{x)> Nhe then no(x) < 1/2 - Mo'e, 
if do(x) < -Mse then uq{x) > 1/2 + Mq' e. 

In view of this last condition, we see that, if is small enough, if do (2;) > 
M2e, then for all < t < /i'^e^l lne|, 

no(x)+e2r(eG, ^) < 1/2 - Mo'e + e^Cg [e^^^^^)! I'^^l/'^ - l] 

< 1/2 + e[ - Mo' + C7ge(M-A'(±£G))/M _ ^Cq] 

< 1/2 + e(-Mo' + 2C6) ^ thanks to (lOnD 

< 1/2 -Cge. 

This inequality and Lemma 15.11 imply ■wf{x,t) < 1/2, where wf is the 
sub-solution defined in (j3.2ip . Consequently, by (|3.25p . 

u''{x,t) < 1/2 if doix) > M2e. 

In the case where do{x) < —M2£, similar arguments lead to u^{x,t) > 1/2. 
This completes the proof of (15. lip . Note that we have proved that, for all 
0<t< /z"ie2|lne|, 

u%x,t)> 1/2 if xGU'i^^XMmA^o), 

u'{x,t)< 1/2 if xen'^0 \^Ahe{ro). 
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Next, since Fj depends on t smoothly, there is a constant C > such 
that, for all t £ [0, /^"^e^l lne|], 

rocAA^,.|i„,i(r,), (5.14) 

and 



(5.15) 



As a consequence of (15.111) and ()5.14p we get 

which completes the proof of Theorem II. 5[ □ 

Proof of Corollary 11.61 In view of Theorem 11.51 and the definition of the 
Hausdorff distance, to prove this corollary we only need to show the reverse 
inclusion, that is 

TtCAfc'em) for 0<t<r, (5.16) 

for some constant C" > 0. To that purpose let C' be a constant satisfying 
C > max(C',C), where C is as in Theorem 11.31 and C as in (|5.15|) . Choose 
t G [0,T], xq € Tt arbitrarily and, n being the Euclidian normal vector 
exterior to Fj at point xq, define a pair of points: 

x^^^ := Xq + e'en and x^^^ := xq — C'en. 

Since C > C and since the curvature of Fj is uniformly bounded as t varies 
over [0, T], we see that, if Eq is sufficiently small, 

Gl^f \AAc;e(FO and x'^'^ e ni'\Mce{rt). 

Therefore, if t G lne|,T], then, by Theorem II. 3 j we have 

u%x^^\t) < 1/2 < u'^{x^ ',t). (5.17) 

On the other hand, if t £ [0,n~'^e'^\ lne|], then from (lEHI) . (I^TB and the 
fact that C > C, we again obtain (pTf]) . Thus (l5l71) holds for ah t £ [0, T]. 
Now, by the mean value theorem, we see that, for each t G [0,T], there 
exists a point x such that 

X G and u'{x,t) = 1/2. 

This implies 5; G Ff . Furthermore we have |a;o — i| < C'e, since x lies on 
the line segment [x^'^^x^^)]. This proves (|5.16|) . □ 
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